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Abstract —Simple RLC circuit models of guitar pickups do 
not account for audible features that characterize the pickup. 
Psycho-acoustic experiments reveal that any acoustically ac¬ 
curate model has to reproduce the first 30 milli-seconds of 
the transient response with extreme precision. The proposed 
model is impractical for simple-minded model reduction or 
brute force numerical simulations yet, by focusing on modeling 
electromagnetic details and exposing a connection to spectral 
graph theory, a framework for finding the transient response to 
sufficient accuracy is exposed. 

Index Terms —Distributed capacitance, spectral graph theory, 
transient response. 

I. Introduction 

The pickup model proposed in this paper has a voltage and 
a current associated with each turn. Although this involves 
thousands of variables, the goal is to devise a model reduction 
strategy which accounts for the audible aspects of the transient 
response which characterize the “tone” [1] of the pickup. The 
“cut-bell” psycho-acoustic experiments of Pierre Schaeffer re¬ 
veal that any acoustically accurate model has to reproduce the 
first 30 milli-seconds of the transient response with extreme 
precision. 1 

The proposed pickup model is not straightforward, but has 
precedent in the analysis of very fast transients in multiwind¬ 
ing power transformers, a field with a history of detailed 
mathematical modeling that goes back a century [3] and where 
the literature on electromagnetics-based state-space analysis of 
transient oscillations, backed by rigorous numerical analysis, 
goes back at least 40 years [4], Problems arising in the context 
of power transformers also shed light on considerations of 
“potting pickups” to suppress microphonics, and of disre¬ 
garding eddy currents on ferromagnetic parts. Unfortunately, 
this analogy does not resolve the issues involved in modeling 
“scatterwound” pickups since the model reduction techniques 
made in the context of power transformers oversimplify the 
pickup problem. For instance, for the purpose of capacitance 
modeling, the turns of transformer winding can be modeled as 
sheet currents [5], [6], [7]. Also, unlike power transformers, 
pickup wire is typically 42 AWG, so skin and proximity 
effects are non-issues well into the MHz range. The key 
common aspect is the neccessity of starting with an accurate 

'Although the brain cannot extract much information based on less than 
three milli-seconds of transient response due to the bandlimited nature 
of the hearing process, it processes transient information quite thoroughly 
based on less than 50 milliseconds of the waveform’s “attack”. Chapter 
12, “Anamorphoses Temporelles I: Timbres et Dynamiques”, of Schaeffer’s 
treatise[2], summarizes this work. There are pointers to earlier observations 
of F. Winckel in this regard, but Schaeffer seems to have provided the key 
experimental confirmation and the readily reproducible methodology. 


multiconductor transmission line model [8], independent of 
questions pertaining to model reduction. 

To fix ideas, consider three questions which one might hope 
to answer in the context of analysing a pickup: 

1) Deriving an accurate transfer function relating the mo¬ 
tion of the string to the output voltage of the pickup; 

2) Finding the input impedance of the pickup. 

3) Find the (short-circuit) internal resonances of the pickup. 

The third problem is a stepping stone to answering the first 
two questions. In this paper we do not attempt to answer the 
first question in full generality. This is because an accurate 
model of pickup would involve, for instance, the eddy-currents 
on the pole pieces of the pickup. Such eddy currents lead to 
“wolf sound”. For a properly adjusted guitar, string heights 
above pole pieces are adjusted to make it imperceptible and 
so we do not need to model the eddy-currents. Similarly, we do 
not dwell on the second problem since the input impedance 
only matters in the context of ensuring that the impedance 
of the guitar and cable, as seen from the follow-on stage of 
amplification, is negligible. Hence, in this paper we focus on 
the role of the of distributed capacitance in the context of the 
third problem. This in turn will lead us to spectral graph theory 
[9], [10], a formalization of the well-studied link between 
nodal analysis, the graph Laplacian and the spectrum of the 
adjacency matrix of a graph. 

II. A Sketch of an Idea 

This paper exposes the connection between the natural 
frequencies of the pickup and the eigenvalues of the adjacency 
matrix. This connection can be sketched as follows. Consider 
the following linear system arising from the discretization of 
multi-conductor transmission lines which are wrapped around 
to form the turns of a winding. It involves square matrices of 
order n, where n is the number of turns: 

(T 2 — Circ)x = 0. (1) 

Here, Circ is a circulant matrix and x are state variables. 
That is, a matrix whose (k+l)-st row is the k-th row shifted 
one entry right and wrapped around. In general any circulant 
matrix is a linear combination of the first n powers of the shift 
matrix S. In other words, Circ is a polynomial in S. Following 
the theory of the DFT, we know S is diagonalized by the 
unitary matrix V, a Vandermonde matrix of roots of unity, 
normalized by n -1 / 2 , and the eigenvalues of S are the n th 
roots of unity. That is, S = V'AgV, where As is the diagonal 
matrix of roots of unity, and f denotes the Hermitian conjugate. 
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Furthermore, all powers of S are also diagonalized by V, as 
are polynomials and analytic functions of S. So, 

Circ = Pi(S) = V t Pi(A s )V (2) 

The square of the propagation matrix f 2 is modeled as a 
polynomial in the adjacency matrix A of a graph associated 
to the manner in which the turns of the pickup winding were 
wound, with the angular frequency w as a parameter: 

r 2 = p 2 (a,w) (3) 

Eq. (3) builds on the relationship between the “Laplacian of 
a graph” where the degree of any node is independent of the 
node index, and the nodal analysis matrix of a network having 
the same underlying graph. In this paper the node degree is 
dependent on the node index at the boundary of the winding 
and an additional diagonal matrix will be introduced to account 
for the discrepancy. Substitute (2) and (3) into (1) to get: 

(P 2 (A,w)-Pi(S))x = 0 or, (4) 

(P 2 (A, w )-V- 1 P 1 (A s )V)x = 0 (5) 

We would expect to find nontrivial solutions by setting the 
determinant to zero to obtain a polynomial equation in io 
whose roots give the natural frequencies. Unfortunately, for 
the model at hand all of the matrices are singular with a 
common null-space and common null vector. Hence the matrix 
is singular for all values of lo. In section six, we will obtain the 
correct equation by using the decomposition in (2) to project 
onto the orthogonal complement of the null-vector, and then 
form the determinantal equation for the natural frequencies. 

III. The Network Model 

Consider a single-coil magnetic pickup on an electric guitar, 
placed in 3-D Euclidean space with coordinates ( x , y, z ) such 
that the z-axis is perpendicular to the fretboard. The pickup 
in question is centered on the z-axis with its center at z=0. 
From a field theoretic vantage, the subtle aspect of the guitar 
pickup model is that each turn of the winding links a nontrivial 
time-varying magnetic field, yet one describes the electric 
field exterior to the conducting coil in terms of an electric 
scalar potential. This is accomplished within the convex hull 
of the pickup by assuming that the time rate of change of 
magnetic flux perpendicular to any plane containing the z- 
axis is negligible. This assumption enables one to define an 
electric scalar potential for the components of the electric field 
lying in any plane containing the z-axis 2 . However, the time- 
varying magnetic flux linking the pickup coil prevents this 
electric scalar potential from describing the circumferential 
components. To describe these components of the electric field 
within the conductor, we integrate the electric field along the 
conductor and evaluate the value of the line integral every time 

2 In principle, the planarity condition in this assumption is only valid for 
a circular coil. A more precise formulation of this assumption would involve 
a singular foliation of the convex hull of the pickup with the all of the 
leaves bordering on the z-axis. A precise formulation of the assumption is 
not required since only the existence of the scalar potential is required to 
formulate the network equations. 


one passes through a given half-plane whose boundary is the 
z-axis. Thus let 

f E-dl = Vi-V 0 (6) 

J a 

where ct comprises the the first i turns of the winding. The 
voltage drop around the f-th turn is then given by V t — Vi- 1 . 
We can compile these voltages into a vector V. Similarly, 
if Ii is the current entering the 7-th winding, these currents 
can be compiled into a vector I. Let n be the number of 
turns. Taking the index modulo n, the equations for the short 
circuited pickup are given by 

(S — \)V = —(R + jwL)7 and, (7) 

(S - I)/ = -jcuCV, (8) 

where I is the identity matrix, S is the shift matrix, and 
the remaining impedance and admittance matrices are as one 
would expect if one took a length of multiconductor transmis¬ 
sion line [8] and connected the end of the fc-th turn to the 
beginning of the (k + l)-st. To preserve the cyclic symmetry, 
we avoided grounding any node and so the capacitance matrix 
is necessarily singular with null vector 1, where 1 is the vector 
all of whose entries are 1. Furthermore, the matrix (S — I) 
has the same null vector. This approach is formalized in the 
method of the indefinite impedance matrix 3 . 

IV. An Analog of the Telegraphist’s Equation via 
Simple Models of R and L Matrices 

If we multiply (7) by jmC we obtain: 

ju>C(S-\)V = -jwC(R+juL)I (9) 

Add and subtract jui (S — l)CV to the l.h.s. of (9) and simplify, 
to get a telegraphist-like equation, plus a commutator term: 

{-ju C(R + ju L) + (S - I) 2 )/ = MCS - SC )V (10) 

Had the capacitance matrix been invertible, and if C and (S —I) 
commuted, we could hope to eliminate V using (8). We will 
see that if we formulate our equations in terms of voltages, the 
current can be eliminated yielding an analog of the telegraphist 
equation, thus avoiding these difficulties. Additional structure 
in the R, L and C matrices leads to useful commutation 
relations. This is easy to verify for the R and L matrices. 
The commutation relations will ensure that the analog of the 
telegraphist’s equations, when formulated in terms of voltages, 
are free of undesirable commutators. 

The coil in a guitar pickup consists of 5,000 to 8,000 
turns of copper wire (typically 42 AWG) with a typical DC 
resistance of 7.5 kfl. The matrix R is modeled as a multiple 
of the identity matrix and is the simplest to model. If /foe is 
the DC resistance of the pickup, and the resistance, ifo of one 

3 A textbook exposition of the indefinite impedance matrix method and its 
relationship to cutset analysis can be found in section 3.7 of Balabanian and 
Bickart’s text [ I]. Note that for a passive network this impedance matrix has 
one zero eigenvalue for each connected component of the network and the 
remaining eigenvalues have positive real pails 
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turn is assumed to be the average, independent of the turn, 
then Ho is Rdc divided by n. Hence, 

R = Rg\ and Trace(R) = Rdc where Rq = —(11) 

n 

The inner turns of the winding are shorter than the outer ones, 
but this is ignored in the analysis. Since R is a multiple of the 
identity matrix, it commutes with any compatible matrix. 

Let M be the mutual inductance between two turns and 
assume this value is independent of the index of the turn. Let 
a be a nonnegative constant expressing the deviation from 
perfect coupling. The inductance matrix model we use is 

L = M(al + (1 §> 1)). (12) 

Here, 1 ® 1 the outer product of the vector 1 with itself 
and again, 1 is the vector all of whose entries are 1. Since 
the capacitance matrix is an admittance matrix arising from 
nodal analysis, and a ground node has not been selected, it is 
symmetric positive semi-definite with a one-dimensional null 
space. This is reflected in the fact that the row sums or column 
sums of any row or column add up to zero. The null vector 
is precisely the vector 1 which occurs in the outer product in 
the inductance matrix. In other words, Cl = 0, so 

C(1 <8> 1) = 0, (13) 

and, by (12) we have 

CL = C(Ma\ + (1 ( 8 ) 1)) = aMC. (14) 

Furthermore, since both L and C are symmetric the same result 
holds for the commuted product, so that 

CL = aMC = LC. (15) 


This argument shows that L commutes with any symmetric 
matrix whose null-space includes the vector 1 and easily 
extends to nonsymmetric matrices where 1 is included in both 
the nullspaces of the matrix and its transpose. One example 
of such a matrix is (S — I), and one can easily verify that 

(S — l)L = aM(S — I) = L(S — I). (16) 


From eqs (7, 8) and commutation relations (15, 16) we obtain: 


(S - \) 2 V 


-(S-l)(R + jwL)I 
-(S - l)(R + juM(a\ + (1 ® 1)))/ 
— (Ro + jujMa)(S — I)/ 
jujC(R 0 + jcuMa)V, or 


((S — I) 2 — jujC(RQ + juMa))V = 0. (17) 


The edges of Gd correspond to neighboring turns, inducing 
a “preferred orientation” on the edges of Gd dictated by the 
order the turns were wound. It follows that the edges of Gd are 
directed from higher to lower node index. Associated to Gd is 
the undirected graph, G u , obtained by ignoring the orientations 
of the edges. We associate the adjacency matrix A to G u by 

... f 0 if i = j or no edge connects nodes i and j , 
(A) ‘J = { 1 otherwise. 

( 18 ) 

If we assume that, locally, the cross-section of the winding 
appears to be a hexagonal closest packing, then any row or 
column corresponding to a non-boundary winding has exactly 
six nonzero off-diagonal entries. The existence of Gd ensures 
that, on average, every node of Gd corresponding to a non- 
bounday turn has three incoming and three outgoing edges. 

From our discussion of the winding’s cross-section, we have 
the following model of the inter-winding capacitiance matrix, 

C = Co(d\ — A+B) and Trace(C) = (nj n td+nbdbav)Co, (19) 

since Trace(A) = 0. Here, C 0 is the capacitance one has 
between two turns of enameled copper wire wound beside each 
other. This parameter is easily deduced by conformal mapping 
once the diameter of the enamel coating is known. The matrix 
B is a diagonal matrix with a relatively small percentage of 
nonzero diagonal entries corresponding to boundary turns. The 
value of the nonzero entry corresponding to a given boundary 
node is the degree of the node minus six. In (19), rii n t and n/, 
are respectively the number of internal and boundary nodes, 
d is the degree of the internal nodes, assumed constant (it is 
six but we leave it as a symbol to suggest that G u is a regular 
graph), and dbav is the average degree of a boundary node. 

Although the incidence data of Gd appears in the statement 
of Kirchhoff’s laws, it is only G u which appears explicitly in 
the nodal equations through the appearance of A. However, in 
the case at hand, the orientation information in the graph Gd is 
implicit in the nodal equations once the preferred ordering is 
introduced. This is because the fc-th edge originates on a node 
whose index j is less or equal to k and terminates on a node 
of index lower than j. In this way, the orientation information 
of Gd can then be reconstructed from the sparsity pattern of 
A, and it is in this manner that winding pattern affects the 
eigenvalue distribution of C through A. 

We have seen that R, L and C are a set of commuting 
matrices and that R, L and (S — I) are another such set. To see 
why C and (S — I) are not, use (19) to reduce the commutator 
on the r.h.s. of (10) to those between A, B and S: 

CS — SC = Co((B — A)S — S(B — A)). (20) 


This is a commutator-free analog of the telegraphist equation. 

V. Winding Induces an Ordering on Turns; 
Adjacency and Distributed Capacitance Matrices 

In order to model the inter-winding capacitance matrix we 
introduce a directed graph, Gd, as follows. Assume that the 
turns are ordered by the order in which they were wound so 
that the i-th index in the matrix equation corresponds to the i- 
th turn wound. This is our definition of a “preferred ordering”. 


Finally, substituting (11,12,14) into (10), reduces it to: 

((S - I) 2 - jcoC{R 0 + juaM))I = jco(CS - SC)V. (21) 

Eq. (20) reveals how formulating our telegraphist equations 
in terms of I leads to unavoidable commutators. And so we 
return to (17) in order to relate the natural frequencies of 
this model to the C matrix via the sparsity structure of the 
adjacency matrix. If we substitute (19) into (17) we get 

((S - I) 2 - juC 0 (R 0 + juMa)(d\ - A + B))V = 0 (22) 
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or 

((S-l) 2 -7 2 (dl-A + B))y = 0; (23) 

7 2 = ju}Co(Ro + jujMa ) (24) 

VI. The tie to spectral graph theory 

Eq. (23) cannot be a conventional eigenvalue problem for 
7 2 since it turns out that the underlying matrices are singular 
for all values of 7 2 . To resolve this problem, we note that the 
vector 1 has appeared in many contexts. Specifically, it is the 
null vector of C = Co(dl — A+B), the null vector of (S —I) and 
the eigenvector of S corresponding to the eigenvalue 1. One 
could toss out a row and column from the system, as when 
one grounds a node in nodal analysis, but a better approach 
is to project the system into the orthogonal complement of 
the vector 1. Recall that the Vandermonde matrix of n-th 
roots of unity diagonalizes any circulant matrix and that the 
eigenvalues of the shift matrix are just the roots of unity. So, 
let Vt denote the Hermitian conjugate of V and let 

An = diag([S2n]£Zo) where fl* = . (25) 

The first term in (23) is a circulant matrix which is diagonal¬ 
ized by the similarity transfomation of section two: 

(S - I) 2 = V t ((A„ - l) 2 )V = V+((A;)V, (26) 

where by (25) we explicitly have 

£ = diag([e^(cos(^) - l)]- 1 . (27) 

Substituting (26) into (23), we can write 

(WA^V — j 2 (dl — A + B))V = 0. (28) 

Although this is still singular for all values of 7 2 , the similarity 
transformation induced by V points the way to projecting into 
the orthogonal complement to the vector 1. Rewrite (28) as 

(A^- 7 2 (dl-V(A-B)V f ))VV=0. (29) 

There is one row and one column which are explicitly zero 
and so we can now explicitly perform the projection. Let 

i '— /r 27 t kj . .27 r/cz. 

A„_i=diag( e " (cos(-) — l)] fc _i- (30) 

n 

That is, with fc ^ 0, n, we restricted ourselves to the nonzero 
eigenvalues of (S — I). Let Vp r be the n x (n — 1) matrix 
consisting of the eigenvectors of (S — I) corresponding to 
nonzero eigenvalues. Then, 

VtV = Wt = l nxn but (31) 

Vp r V^ r = l(„— 1 ) X („— 1 ) and V}, r V Pr . = Proj _L 1 (32) 

The determinantal equation is now expressed as: 

det(A^Ti — 7 2 (c(l(n-i)x(n- 1 ) — Vp r (A — B)Vp r )) = 0 (33) 

This is the definitive equation for subsequent developments. 
However, given the definition of the capacitance matrix, this 
equation could have been written in the more intuitive form 

det(vCTi + ^-Mpr CV^.) = 0. (34) 

Go 


Hence, (33) reduces the problem to a modified eigenvalue 
problem for a (diagonal) perturbation of the (vanishing diag¬ 
onal) adjacency matrix. By (24) any information about the 
spectrum translates into information about the collection of 
natural frquencies since they are related to the eigenvalues 
through a quadratic equation. The adjacency matrix eigenvalue 
problem is the focus of “spectral graph theory”, and there are 
many ways to exploit this connection. This provides precise 
estimates for the transient response and a framework for model 
reduction in the presence of eigenvalue clusters. 

VII. Summary and Conclusions 

In this paper we shed light on elusive aspects of electro¬ 
magnetic modeling of electric guitar pickups. Specifically, the 
ability of an expert to infer the manner in which the pickup 
was wound by listening. Psycho-acoustic experiments reveal 
that acoustically accurate models have to reproduce the first 30 
milliseconds of the transient response with extreme precision 
since the brain makes inferences about the pickup’s tone within 
this timeframe. Conventional wisdom dictates that one could 
take the natural frequencies which are smallest in absolute 
value and perform a model reduction. This however is tricky 
when the natural fequencies cluster. The problem is difficult 
because the winding pattern of the pickup affects the clustering 
of natural frequencies in a subtle manner which is audibly 
discernable as the tone of the pickup, yet gets lost in a simple- 
minded model reduction scheme. 

By exposing a connection to spectral graph theory, 
(23,24,33) yield a framework for analyzing the “attack” of the 
transient response to sufficient detail to reveal the role of the 
ordering of the turns in the winding. This model is insensitive 
to how the pickup was wound with the exception of the inter¬ 
winding capacitance matrix which describes the capacitance 
between the turns of the winding. 
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